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ABSTRACT 

We present the results of extinction measurements toward the main ejecta 
shell of the Cassiopeia A supernova (SN) remnant using the flux ratios between 
the two near-infrared (NIR) [Fe II] lines at 1.26 and 1.64 pm. We find a clear 
correlation between the NIR extinction {E{J — H)) and the radial velocity of 
ejecta knots, showing that redshifted knots are systematically more obscured 
than blueshifted ones. This internal “self-extinction” strongly indicates that a 
large amount of SN dust resides inside and around the main ejecta shell. At 
one location in the southern part of the shell, we measure E{J — H) by the 
SN dust of 0.23 ± 0.05 mag. By analyzing the spectral energy distribution of 
thermal dust emission at that location, we show that there are warm (~ 100 K) 
and cool (~ 40 K) SN dust components and that the latter is responsible for 
the observed E{J — H). We investigate the possible grain species and size of 
each component and find that the warm SN dust needs to be silicate grains such 
as MgSiOs, Mg 2 Si 04 , and Si02, whereas the cool dust could be either small 
(< 0.01 pm) Fe or large (> 0.1 pm) Si grains. We suggest that the warm and 
cool dust components in Cassiopeia A represent grain species produced in diffuse 
SN ejecta and in dense ejecta clumps, respectively. 

Subject headings: dust, extinction — infrared: ISM — ISM: individual objects 
(Cassiopeia A) — ISM: supernova remnants 


1. Introduction 

For the past decade, evidence for copious amounts of dust (> 10® Mq) in high-redshift 
galaxies and quasars with z > 5 has been reported in many far-infrared (FIR) and sub- 
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millimeter/millimeter studies fe.g.. lWang et al.ll201ll:iLeipski et al.ll2013l: ICalura et al.ll20l4 
and references therein), implying that dust should be produced within a very short timescale, 
less than 1 Gyr after the Big Bang. The most promising sites of dust formation in such an 
early universe are believed to be core-collapse supernovae (CCSNe), since their progenitors 
have very short lifetimes (~ 10® years) comp ared to asymptotic giant branch stars, whic h 
are major sources of dust in the local universe (IMorgan fc Edmnndsll2003l: lMarchenkoll2006l) . 


To understand the dust formation process in expanding; supernova (SN) ejecta, many _ 

theor etical studies have been conducted to date (INozawa et ahll2003l . l2008l. 120101: ISarangi fc Cherchnefl 
2014J) . According to these theoretical studies, dust can form in an expanding He core with 
grain species determined by the elemental composition at the formation site. If an SN ejecta 
maintains its layered structure during expansion, a variety of grain species can form, e.g., 
carbon grains in the helium layer, silicate and oxide grains in the oxygen-rich layer, and 
Si/Fe grains in the innermost layer, the amounts of which depend on the progenitor mass. If 
the SN ejecta is mixed, however, it is mostly silicate and oxide grains that can form. The size 
of the grains depends on the gas density at the condensation time. In an SN IIP, which has 
a thick hydrogen e nvelope, the expansi on is slow, so large grains with an average radius of 
~ 0.1 yum can form flNozawa et al.ll2003l) . In an SN Ilb with a low-mass hydrogen envelope as 
well as envelope-stripped SNe such as Type Ib/c, on the other hand, th e SN ejecta expands 


rapid ly, so only small grains with a size less than 0.01 /mi can form flNozawa et al.l 12008 


f 


affected bv densitv enhancements due to clumping fSarangi & Cherchnef 

2014 

). Theoreti- 

cal studies predict that one CCSN can produce 0.1-1.0 Mp. of SN dust (h 

fozasa e 

al. 1991: 

Todini & Ferrara 2001: Nozawa et al. 2003. 2008. 2010: Bianchi & Schneider 2007 

), but the 


dust masses obtained from observations of nearby CCSNe are much smaller than this, e.g., 


< 10 ^ Mq f Meikle et ahl 120071: IXotak et al.ll2009l: ISakon et al.ll2009l: ISzalai et al.ll2011 h A 


remnants (SNRs), e.g., Cassiopeia A (Barlow et a 

2010: Sibthorpe et al. 2010l: Arendt et al. 

20141. Crab (Gomez et al. 2012: Owen & Barlow 

2 OI 5 I 1 . and 1987A (Matsnnra et al. 2011: 

Indebetouw et al. 2014ll. but the phvsical and chemical properties of the dust remain to be 


explored (see below). 


Cassiopeia A (Cas A), which is a remnant of an SN lib explosion flKrause et al.ll2008l) . 
is an ideal target for examining the properties of SN dust. It is quite young (~ 330 years; 
Fesen et al.ll2006l) . so the SN ejecta and freshly formed SN dust have not been signihcantly 
mixe d with the circum stellar/interstellar medium (CSM/ISM), and it is relatively close (3.4 
kpc; Reed et al. lOOsI) . so the physical and chemical structures can be resolved. Therefore, 
Cas A has been a major target of infrared space missions. The hrst direct evidence for SN 
dust in Cas A came from the Infrared Space Observatory [ISO), which performed mid-infrared 


























































































































3 


(MIR: 2.4-45 ^ni) spect r oscopic observations toward the bright ejecta shell flLaeage et ah 
19961: lArendt et ahlll999l: iDouvion et al.l 1200111 . The continuum spect ra had a strong bump 
peaking at 21 /rm together with a relatively weak bump at 9.5 /rm, and iDouvion et ahl ( 20011) 
showed that the spectra can be htted by two dust components at different temperatures, one 
at 90 K and the other at 350 K, with pyroxene (MgSiOs), quartz (Si 02 ), and aluminum 
oxide (AI 2 O 3 ) as major components. The composition and distribution of this “warm” SN 
dust heated by a reverse shoc k have been studied in detail using the spectral mapp ing data 
of the Spitzer Space Telescope ( Ennis et ah 2006 : Rho et ah 2008 : Arendt et ah 2014 ). These 
studies showed that the dust emission exhibits distinct spectral characteristics depending on 
the SN material with which it is associated. The most prominent dust emission is that with 
strong 9 and 21 fim bumps, which is associated with SN ejecta having strong Ar emission 
lines. The spectral features agree with the ISO spectra and can be reproduced by Mg 
protosilicate/MgSiOs, Mgo. 7 Si 02 . 7 , or Si 02 . The dust emission associated with SN ejecta 
having strong Ne lines is smooth without any silicate features and can be reproduced by 
AI 2 O 3 dust or C glass. There is yet another type of smoo th dust emission as sociated with 
X-ray Fe emis s ion, b ut it is probably from swept-up CSM (lArendt et ahl 120141) . Meanwhile, 
Nozawa et al.l (120101) showed that the observed spectral energy distribution (SED) of dust 
emission can be fitted by a physical model in which the SNR is expanding into a dense CSM, 
and the MIR emission is mostly from silicate grains (e.g., MgSi 03 , Mg 2 Si 04 , and Si 02 ) and 
MgO. The estimated total mass of warm dust ranges from 0.008 to 0.054 Mq. 

A larger amount of SN dust, however, appears to reside in th e interior of Cas A, where 
the SN ejecta is freely expanding. The first report was made bv iDnnne et ahl (120031) . who 
claimed the detection of 2-4 Mq of cold (15-20 K) SN dust associated with Cas A from an 
850 fim observation. It was found, however , that at least some of the 850 /im emission is from 


a fore groun d material (iKrause et ahl |2004 iDunne et al.l 120091) , and instead ISibthorpe et ah 


(I 2 OIOI) and iBarlow et al.l (120101) . using the balloon-borne sub-millimeter telescope BLAST 
and infrared space telescopes AKARI and Herschel, showed that there is faint FIR and 
sub-millimeter emission from Cas A that can be fitted by dust at 35 K with a total mass 
of ~ 0.07 Mq. lArendt et al.l (120141) also identified the FIR emission from this “cool” dust 
associated with unshocked SN ejecta (their “[Si II]” dust group) using the HerschelFIR data 
and derived an upper limit of 0.1 Mq for the dust mass. However, the composition of the 
cool dust has been undetermined so far because the absorption cross sections of dust grains 
in the FIR are mostly smooth without prominent spectral features. 

In this paper, we show that one can explore the characteristics of the cool SN dust in 
the interior of Cas A using extinction in the near-infrared (NIR). Recently, we performed 
NIR spectroscopy toward the main ejecta shell of Cas A, where we obt ained the spectra l 
and kinematical properties of 63 knots bright in [Fe II] emission lines (iKoo et ahl 120131) . 
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Most of these knots are SN material, and the extinctions toward these knots show clear 
evidence for the SN dust, which can be combined with thermal infrared dust emission to 
infer the dust composition. This paper is organized as follows. In Section [21 we outline 
our NIR spectroscopic observations and data reduction and then briefly describe the public 
MIR and FIR imaging data that we used in the paper. In Section [3l we derive the NIR 
extinction toward the knots using [Fe II] line ratios and show that there is a correlation 
between the extinction and the line-of-sight velocity, which implies extinction by SN dust. 
Then, in Section 01 by comparing the NIR extinction and the thermal dust infrared SED, we 
constrain the possible compositions and sizes of the SN dust. We discuss the implication of 
our result for dust production in Cas A in Section |5] and summarize our paper in Section |6l 


2. Observations and Data Reduction 


2.1. Near-infrared Spectroscopy 


Eight long-slit NIR spectra of the main SN ejecta shell were obtained on 2008 June 29 
and August 8 using TripleSpec mounted on the Palomar Hale 5 m telescope; this spectro¬ 
graph provides a broadband spectrum simultaneously covering 0.94-2.46 /rm with a moder¬ 
ate resolving power of ~ 2700. We also observed the spectra of an AOV standard star (HD 
223386) for flux calibration just before and after the target observations. To subtract sky 
OH airglow emission as well as background continuum, the observations were made with 
either a consecutive ABBA pattern or On-Off mode depending on the complexity of the 
held. Figure [T](a) shows the position of the slits on the continuum-subtracted [Fe II] 1.64 pm 
narrow-band image obtained in 2008 August. 


We followed a general reduction procedure to make a two-dimensional dispersed image. 
Numerous sky OH airglow emission lines were used for wavelength calibration, and then 
we corrected it to the heliocentric reference frame. The absolute photometric calibration 
was made by comparing the observed spectrum of the standard star with the Kurucz model 
spectrum f)Kuruczll20031) . Even though the absolute hux uncertainty goes up to 30% depend¬ 
ing on the centering accuracy of the standard star observation, the relative huxes are quite 
robust. 


We found a total of 63 infrared knots using a clump-hnding algorithm (IWilliams et ah 


I 994 J ) and performed a single Gaussian £t for all the detected emission lines to derive their 
huxes and line-of-sight velocities. For detailed descriptions of the data reduction a nd analysis 
as wel l as the full list of parameters of the detected emission lines, please refer to iKoo et ah 
(12013!) and Y.-H. Lee et ah (2015, in preparation). 
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2.2. Mid- and Far-infrared Imaging Data 


For an SED analysis of SN dust, we fully exploit the previously published archival 
data providing the highest spatial resolution and sensitivity in the MIR and FIR wavebands 
covering 10-500 /ini, where the dust emission is dominant. 

For short wavebands, we us e the W3 channel i mage (Aiso ~ 12 pm) of the Wide-field 
Infrared Survey Explorer ( WISE: IWright et al.ll2010h with an angular resolution of 6''5. We 
hrst retrieved all 36 single-exposure images (“Level-lb”) from the online servic^ and con¬ 
verted the pixel unit to MJy sr“^. After subtracting the predicted brightness of zod iacal 
light calculated from the interplanetary dust model from each frame flKelsall et al.lllOOSh . we 
co-added the individual frames using Montage softwar^ to obtain the hnal image. Instead 
of the WISE W4: channel image (Ajgo ~ 22 /im) , we use the 24 pm band image of the Multi¬ 
band Imaging Photometer for Smtzer (MIPS: iRieke et al.ll2004l) obtained in 2007 January 
(AORKEY: 17657088, 17657344 )i, which provides better spatial resolution (6") than WISE 
W4. The data used here are the Post-BCD image (“Level 2”), which is automatically pro¬ 
cessed by the MIPS software pipeline (version; S18.12.0). As in the WISE data, zodiacal 
light is one of the signihcant sources of background radiation in this waveband. We therefore 
subtracted the bac kground brightness from the mosaic image using the model calculation 
( Kelsall et al. 1998 1. 


Recent FIR images of t he remnant were ob served with PACS and SPIRE on board 
Herschel Space Observatory fjPilbratt et al.l 1201011 in 2009 September and December. The 
two instruments have six photometric bands (70, 100, and 160 pm for PACS, and 250, 350, 
and 500 pm for SPIRE), and they provide spatial resolutions of 572, 777, 12", 18", 25", 
and 37", respectively. The fully calibrated data sets were downloaded from the Herschel 
Science Archiv^, and we use the data product levels of 2.5 and 2.0 for PACS and SPIRE, 
respectively. All imaging data have well-established zero-level brightness, whereas the PACS 
photometer images at level 2.5 do not (PACS Observer’s Manual - version 2.5.1). To obtain 
the true zero point of the images, we calibrated the zero level by comparison with refer¬ 
ence images whose zero point is well established. For the reference photometric calibrator, 
we used 60 and 100 ^m images fr om Improved Reprocessing of the IRAS Survey (IRIS: 
Miville-Deschenes fc Lagachel 120051) together with DIRBE 140 and 240 pm and Planck 350 


^http://irsa.ipac.caltech.edu/applications/wise/ 

^http://montage.ipac.caltech.edu/ 

^http://sha.ipac.caltech.edu/applications/Spitzer/SHA/ 
^http://www.cosmos.esa.int/web/herschel/science-archive/ 
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and 550 /im maps. In Cas A, synchrotron radiation is a dominant source in the FIR wave- 
band. This no nthermal component is subtracted by using the VLA radio image at 6 cm 
flDeIanevll2004l) . The expected brightness at each isophotal wavelength is extrapolated from 
the power law 8^ = 0^° (erg cm“^ s“^ Hz“^). In this paper, we adopt a = —0.682 and 
C = 2.50 X 10® erg cm“^ s“^ Hz“^, which a re derived from the total brightnes s of the remnant 
at 30, 44, 70, and 100 GHz using Planck f Planck Collaboration et al. 2011) archival data@. 


3. NIR Extinction by SN Dust 
3.1. NIR Extinction Measurement 

We derive the selective extinction or “color excess” between 1.26 and 1.64 /mi, E{J—H), 
toward the infrared knots from 


E{J — H) = Ai .26 — Ai ,64 = 1.086 In 


[F’i. 26/ -^1 .eijint 

[Fi.2,IF, 

.eijobs 


( 1 ) 


where Ai ,26 and Ai .64 are the total extinctions at 1.26 and 1.64 /im, and [Ei. 26 /-^i. 64 ]obs and 
[-^i. 26 /-^i. 64 ]int are the observed and intrinsic flux ratios of the [Fe II] 1.26 and 1.64 /im lines, 
respectively. These two strong [Fe II] lines share the same upper state (a'^D 7 / 2 ), so their 
intrinsic ratio is hxed by their Einstein A coefficients and wavelengths, i.e., [Ei. 26 /-^i. 64 ]int = 
(Afci4.26/F26)/(Afcj_i,64/1.64). The ratio has been derived from both theoreti cal calculations 
and observations, and it ranges from 0.94 to 1.49 in previous studies (e.g., iGiannini et al. 
20151: IKoo fc Leell2015l. and references therein). This large uncertainty of the intrinsic flux 
ratio hampers accurate measurement of the absolute extinction toward the remnant (this 
issue will be addressed in Section 


It is, however, worth noting that the uncertainty 
does not affect the difference in E{J — H) (Equation d]) among the knots, which gives the 
relative extinction to them. In this pap er, we adopt the intrinsic lin e flu x ratio of 1.36, 
which is the theoretical value proposed by iNussbaumer fc StorevI fjl988|) and iDeb fc Hibbert 

fcoioh . 


3.2. Total Extinction to Cas A 


In Figure [I][b), we plot E{J — H) for the 63 infrared knots as a function of position 
angle. The colors and symbols represent the characteristic groups classihed by iKoo et al. 


“http://irsa.ipac.caltech.edu/applications/planck/ 
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(120131) : helium-rich knots in green squares, sulfur-rich knots in red circles, and iron-rich 
knots in blue diamonds. In short, helium-rich knots are slowly moving (< 100 km s~^) knots 
with strong He I lines, the properties of which match those of “qua si-stationary flocculi' 


(QSFs), the ci rcumstellar material swept up by an SN blast wave (Ivan den BerghI Il971 
Lee et ahl 12014) . Sulfur-rich and iron-rich knots are fast-moving (> 100 km s“^) knots with 


strong [S II] and strong [Fe II] lines, respectively. They are SN ejecta material synthesized in 
different laye rs of the SN and correspond to fast-mqying k nots (FMKs) in previous optical 
studies fe.g.. lvan den BerghI Il97ll: iHammell fc FesenI 120081) . For more detailed ex planations 


on th e spectroscopic and kinematic properties of these knots, please refer to iKoo et al. 


r 

Ik 


(120131) . For compari son, we also plot 7 71J — H) for the optical and infrared knots measured 
in previous studies. Hurford fc Fesen ( 19961) observed optical spectra of hve FMKs and two 
QSFs located in the northern bright rim and measured E{B — V) using ratios of the [S II] 
and the Balmer lines, respectively. We converted it to E[J — H) (= 0.30 x E (B — V)) 
assuming the general interstellar dust composition with Ry of 3.1 (iDraind 120031) . and the 
open red circl e s and green squares in the hgure represent the FMKs and QSFs, respectively. 
Eriksen et al.l (120091) performed NIR spectroscopy for four slit positions and derived the line 
ratio of [Fe II] 1.26-1.64 /im for 19 infrared knots without any classihcation of their origin. 
We obtained E{J — H) using Equation [T] and marked the results by black X marks. 

The derived E{B — V) varies considerably over the remnant. It appears that the extinc¬ 
tion toward the west is systematically larger than that toward the east, though there is a large 
scatt er. This systematic variation is well known f rom previous radio and X-ray observations 
(e.g., Keohane et al. 19961: Revnoso fc Goss 2002 : Hwang fc Lamin j2012|j. For example, the 
green contours in Figure[T](a) show the column density (A^n) map of lHwang fc Lam in e 



obtained from an analysis of the Chandra X-ray data. The column density is large ((2-3) 
xlO^^ cm“^) toward the west and southwest directions, and molecular line studies showed 
that it is due to molecular clouds located in the Perseus Spiral Arm. 

Figure [11(b) shows that the extinction varies considerably within a slit. For example, 
along Slit 4, E{J — H) varies from 0.6 to 1.5 mag even though the emission is from a thin 
hlament. This variation is signihcantly larger than what we would expect from the variation 
in the foreground extinction. For example, the solid line in the hgure represents an expected 
E{J — H) variation along the ejecta shell derived from the Ah niap in Figure [T](a). The 
ejecta shell is determined in the continuum-subtracted [Fe II] image where the brightness 
is greater than 3cr of the background rms noise, and the mean Ah is obtained at every 10° 
in position angle. The column density is convert ed to E{J — H) using Ny{/E{B — V) = 
5.8 X 10^^ atoms cm“^ mag“^ ( Bohlin et al. 1978 1 and E{J — H) = 0.30 x E{B — V) = 


Ah/1. 9x 10^^ cm The dashed-dotted lines represent the maximum and minimum E{J—H) 
at each position angle. It is obvious that, in several slits, the variation oi E{J — H) within a 

























































slit is much larger than the range allowed from the Nu map. As we show in the next section, 
this large variation is due to extinction within the SN ejecta. 


3.3. Self-extinction within Cas A 


An interesting correlation is found when we plot E{J — H) as a function of radial 
velocity, which is shown in Figure [2]^ a). It is clear that the redshifted knots are generally 
more obscured than the blueshifted knots. Because we are observing an expanding shell, this 
suggests that the knots on the far side experience more extinction than those on the front 
side, or that there could be extinction originating within the remnant. To conhrm this “self¬ 
extinction,” it is necessary to remove the spatially varying foreground interstellar extinction. 
We subtract the extinction derived from the X-ray observation (Figured]), and the result is 
shown in Figure |2](b). The correlation becomes slightly weaker, but it is still obvious that the 
redshifted knots are more heavily obscured than the blueshifted knots. Therefore, our result 
implies that there is extinction, which might be due to newly formed SN dust within the 
remnant. The large scatter in extinction for the same velocity knots might reflect differences 
among different sight lines, i.e., non-uniform distribution of SN dust. 


Note in Figure[2](b) that the minimum E{J—H) ~ —0.3, which means that the E{J—H) 
derived from the Ah map is higher than the E{J — H) derived from [Fe II] line ratios. We can 
think of two possibilities for the negative E{J — H). First, it could be due to the uncertainty 
of the theoretical [Fe II] line ratio. To compensate for the E{J—H) of —0.3 mag, however, the 
theoretical flux ratio of [Fe II] 1.26-1.64 pm should be more than 1.8, which is considerably 
higher than the numerical value obtained in previous studies (0.94 < [Ei^ 26 /Fi.m]int < 1.49). 
Second, the Ah derived from X-ray data could be overestimated. Indeed , the mean column 
densi ty from the X-ray studies is ~ 1.5 x 10^^ cm“^ flWillingale et al.ll2002l:lHwang Sz Laming 


1985 


2012 ) , which is higher thaii those of radio observations, e.g., ~ 1.1 X 10^^ cr n ^ flTroland et ah 


1I2OI2I) 


Keohane et al.Nl996[) . We also note that the colum n densitie s of iHwang Sz Laming 


are higher by (0.4-0.7) xlO 


22 —2 
cm ^ 


than those of iLee et al.l (120141) . who carefully 
analyzed the X-ray spectra aronnd the onter SNR shock. The difference of 0.3 mag in 
E{J — H), which corresponds to Ah ~ 6 x 10^^ cm~^ , may therefore have res ulted from 
the difficulty in backgronnd removal in the analysis of iHwang fc Laming] (120121) . Another 
complication is that the X-ray-absorbing colnmn also inclndes self-extinction, so the X-ray- 
based extinction overestimates the foregronnd extinction. Thns it appears difficult to obtain 
an accurate map of the foreground extinction from the X-ray data alone, but this uncertainty 
is not likely to erase the systematic correlation in Figure |2][b). 


The non-uniform spatial distribntion of SN dust and the uncertainty in the foregronnd 
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extinction hamper the analysis of the correlation in Fignre |2t^b). However, in one slit (Slit 
4), we detected convincing evidence for self-extinction withont those sight-line-dependent 
complications. Fignres ini(a) and (b) show that the velocity strnctnre of the [Fe II]-hne- 
emitting gas in Slit 4 is “arc”-like, snggesting that the slit crosses a portion of an expanding 
shell. This is consistent with the resnlt of pre vious stndies that the main ejecta shell is a 
thin shell expanding at 4000-5000 km s~^ fe.g.. IPeLanev et al.ll2010l: llsensee et al.ll2012l) . In 
Slit 4, therefore, we have both blue and redshifted Fe knots along a given sight line, and the 
uncertainty due to different sight lines disappears. Now Figure[21(c) shows the E{J — H) map 
obtained for pixels with an [Fe II] 1.64 /rm brightness greater than 3cr above the background 
noise. A systematic increase in the extinction toward the redshifted knots is apparent. We 
computed an intensity-weighted E{J — H) as a function of radial velocity, which is shown by 
the black solid line in Figure 131(d). We also divided the iron hlament into three subregions 
along the slit (A-C in Figure [31(c)), and their mean E{J — H) values are over-plotted on the 
hgure with red, green, and blue solid lines. The total E{J — H) within the ejecta shell from 
Figure [31(d) is 0.23 ± 0.05 mag with a systematic gradient of 0.13 ± 0.02 mag per 1000 km 
s“h This should be mostly, if not entirely, due to SN dust because the blue and redshifted 
ejecta materials in the slit are essentially along the same sight line. There are small-scale 
variations in E{J — H) in Figure [31(d). They could be due to either incomplete subtraction 
of OH airglow emissions (e.g., gray-hatched areas in the hgure) or very small dust clumps 
within the slit, or both. 


The obtained E{J — H) can be converted to the dust column density (Ed) with an 
appropriate Akjh = Ki .26 — fi^i.64 
coefficients at 1.26 and 1.64 /rm, respectively: 


cm^ g ^), where Ki ,26 and are the mass absorption 


E{J — H) = 1.086 / Akjh Pd dl ~ 0.11 


Ak, 


JH 


10^ cm^ g ^ 


10“"^ g cm“2 


( 2 ) 


This equation yields the column density of a given dust species required to explain the 
observed NIR self-ext inction. For example, if it is Fe dust, which has Akjh ~ 10^ cm^ g 


-1 


fISemenov et ahl[20031. see Figure [7D, the required dust column is ~ 10 ^ g cm 


4. FIR Emission and SN Dust Composition 
4.1. FIR Emission from SN Dust 

What dust species would produce the E{J — H) of 0.2 mag detected in Slit 4? In this 
section, we constrain the dust species by investigating the FIR emission properties of SN 
dust. 














- 10 - 


Figure m^a) shows the one-dimensional (ID) brightness prohles at 12, 24, 70, and 100 
/rm along the slit length in Slit 4. Note that these data have similar angular resolutions of 
5''2-7'/7 fSection I2.2p . The prohles are background-subtracted using the mean brightnesses 
of an area just outside the ejecta shell (red square in Figure 111(b)) and are normalized by 
their peak brightnesses. All four prohles have maxima at 1'.67 where the iron hlament is 
located (Figure |3]). All emission drops to zero inside of the ejecta shell except for the 100 
fim emission, which remains constant at ~ 30% of its peak brightness. This excess emission 
at 100 fim must be mostly from the un shocked cool SN dust detected previously in the FIR 
( Barlow et ah 2010 : Arendt et ah 2014). 


We examine the SEDs at two positions (see Figure 111(a)): (1) the peak brightness position 
at 1(67 with ~ 6" aperture (hereafter position A) and (2) the inner 100 pm excess region 
between 1(2 and 1(4 (hereafter position B). For comparison, we also derive an interstellar 
dust SED just outside of the shell (red square in Figure 111(b)) where the radiation from 
the remnant is almost negligible. To derive accurate parameters for the SN dust, we hrst 
have to subtract the contributions of line emissions from the measured brightnesses. For 
this purpose, we retrieved the archival MIR and FIR spectra from Spitzer IrM and ISO 
LW^ observations , respectively. We produced an IRS spectral data cube using CUBISM 
flSmith et al.ll2007l) and extracted ID spectra within a 6 " aperture at positions A and B as 
well as the background position. The strong emission lines falling in the WISE W3 and 
MIPS 24 pm bands are the [Ar II] 8.99 pm, [S IV] 10.5 pm, [Ne II] 12.8 pm, [Ne III] 15.6 
pm, and [O IV]-|-[Fe II] 25.9 pm lines, and their contributions to the fluxes in the two bands 
have been estimated as 8.9% and 4.8% for position A, 4.6% and 16% for position B, and 
2.0% and 1.8% for the background position, respectively. There are no corresponding FIR 
spectroscopic data, but an ISO LWS I D spectrum (43-190 um) f rom an area close (~ 50") 
to Slit 4 is available (#3 spectrum of iDocenko fc SunvaevI (120101) 1. According to this ISO 
spectrum, the strong emission lines falling in the PACS 70 and 100 pm bands are the [O I] 
63.2 and [O III] 88.4 pm lines, respectively. Estimation of the contribution of these lines to 
the observed surf ace brightne sses at our positions is tricky because the ISO LWS beam is 
large, i.e., ~ 80" flLlovdl 120031) . and covers both the Cas A ejecta shell and the background 
area. We estimated the line contribution by assuming that both emission lines come from the 
ejecta shell, not from the interstellar material. This assumption is plausible because the two 
lines are detected at several positions along the ejecta s hell to show similar broad prohles, 
whereas there is little emission outside of the remnant (IDocenko fc SunvaevI l2010l ). Under 


this assumption, the two lines contribute only to the brightness of A, which is estimated as 


®http://sha.ipac.caltech.edu/applications/Spitzer/SHA/ 
^http://iso.esac.esa.int/ida/ 
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3.8% and 4.5% in the PACS 70 and 100 /im bands, respectively. In the PACS 160 /im band, 
there is the [O I] 146 /im line, but its contribution is almost negligible. The observed and 
line-subtracted dust continuum brightnesses at the three positions are listed in Table [U 

Figure |5] shows the line-emission-subtracted SEDs at positions A and B in red and blue 
solid lines, respectively. Both SEDs have broad thermal emission bumps peaking at 160 
/im, which arise from the line-of-sight interstellar dust component at a temperature of ~ 20 
K. This FIR SED matches well the SED of the background area outside the shell (black 
solid line), which represents the interstellar dust SED. The background-subtracted SEDs, 
which presumably represent “pure” SN dust SEDs, are shown in red and blue dotted lines 
in Figure |5l The SED at A now shows maximum brightness at 24 /mi, and the brightness 
decreases toward longer wavelengths. On the other hand, that at B peaks at 100 /im, 
which implies that the dust temperature is much lower than that at A. These two SEDs 
should represent the SEDs of shocked (-|-unshocked) and unshocked dust species, respectively. 
We will analyze the SED at position A, where we obtained self-extinction of E{J — H) = 
0.23 ±0.05. 


Before doing a detailed analysis in the next section, we hrst show that the shocked 
warm dust component emitting the MIR emission cannot explain the NIR color excess 
that we detected. Previous observations showed that the warm dust is silicate grains (e.g.. 


Arendt et al.l 120141) . and we assume that it is MgSiO.s grains w ith a size of 0.001-0.1 /im, 

g“^ flDorschner et al.lll995l) . (As we will show in the 


which has 10 ^ 5 -, Ak 


JH 


< 10^ cm^ 


next section, assuming other silicate dust species will yield the same conclusion.) The SED 
at A has excess emission in the FIR, and we perform a two-temperature modihed black- 
body ht assu ming t hat the cool dust has t he optical properties of general interstellar dust 
(iDrainel 120031) . as in ISibthorpe et al.l (120101) . The surface brightness of each component at 
temperature T is calculated from 


Su = Kabs,u By{T) Ed 


(3) 


where is the mass a bsorption coefficient, and By{T) is the Planck function. We used 

MPFIT flMarkwardtll2009l) . which is a nonlinear least-squares htting routine within IDE, to 
derive the values of T and Ed at which is minimum as well as their formal la uncertainties. 
The SED is well htted with warm dust at T = 115 ± 2 K plus cool dust at T = 40 ± 6 K 
(Figure [5]). From this two-component htting, we found that almost all of the brightnesses 
at 12 and 24 /rm arise from the warm dust component, whereas more than 80% of the 
brightness in the FIR longer than 100 /rm is from the cool dust component. The warm and 
cool dust components yield Ed = 3.2(±0.3) x 10“® g cm“^ and 1.1(±0.8) x 10“® g cm“^, 
respectively. From Equation[2l therefore, the corresponding E{J — H) of the warm dust with 
a size of 0.001-0.1 /mi is 10“® — 10“^ mag, and that of the cool dust is ~ 0.003 mag. The 
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column density of warm dust required to explain the MIR emission is very small, and the 
resulting NIR extinction is several orders of magnitude smaller than the value we derived 
(0.23 ±0.05 mag). Therefore, it must be the cool dust component that is responsible for the 
NIR extinction. (Note that the grain species of the warm dust component and therefore its 
contribution to the NIR extinction are rather well constrained.) However, the above result 
suggests that the cool dust component cannot have the optical properties of the general 
interstellar dust because then the NIR extinction is again two orders of magnitude less than 
the value we derived. The grain species of the cool dust should have large Akjh and/or 
small FIR opacity (^abs,!/) compared to that of general interstellar dust. 


4.2. SN Dust Composition 
4-2.1. Dust Species and Their Optical Properties 


For the composition of S N dust, we consider th e grain species predicted by pr evious 
theoretical models for CCSNe ( Todini fc Ferrara 2001 : Nozawa et ah 2003 . 2008 . 2010 ). The 
grain species considered in this paper are summarized in Table 2 with references for their 
optical constants {n and k). Note that carbon grains are formed in the He layer, silicates 
(MgSiOa, Mg2Si04, and Si02) and oxides (AI2O3 and MgO) in the 0-rich layer, and other 
heavy-element (Si, Fe, FeS) grains in the Si-S-Fe layer. In addition, various types of oxide 
grains (SiO^, MgSiO.s. Mg2Si0 4, AI2O3, and Fe304) are produced in the mixed SN models. 
According to iRho et al.l (120081) . the MIR spectra of Gas A are well htted by the combination 
of these grain opacity curves, which implies that these are the major dust species produced 
by the SN explosion. We cal culated the absorpt i on an d extinction coefficients (Kabs.jy and 
^ext.i/) using the Mie theory (iBohren fc Huffman! 119831) assuming a spherical grain with a 
radius of 0.001, 0.01, or 0.1 pm. 

Figure E] shows Akjh and Kabs,i/ at 24 and 70 pm for the grain species with different sizes. 
As expected, the Kabs,u values of dust species at 24 and 70 pm do not depend on the grain 
size except for Fe dust. For Fe dust, which is a metallic dust, not only the electric dipole but 
also the magnetic dipole is important for grains larger tha n 0.006 um, the effi ciency (Qabs) 
of which is proportional to the volume and not the size flTanabe et al.Nl984i) . Therefore, 
the FIR opacity (Kabs,^ = (7ra^Qabs,i^)/(4/37ra^p), where Qa.hs,u is the absorption efficiency) 
of Fe dust is proportional to the square of the size for large grains. Note that small Fe dust 
(< 0.01 pm) has an FIR dust opacity much smaller than those of the other dust grains, 
e.g., < 0.2 cm^ g“^ versus 2-100 cm^ g“^ at 70 pm. Anju in general does not depend on 
the grain size for a < 0.01 pm and increases with a for larger grains up to a > 0.1 pm 
(Figure [7j). This is the result of standard grain extinction properties: for small grains, i.e.. 
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X = 2'7ia/X <C 1, Qext(= Qahs + Qsca) IS dominated by Qahs, which is proportional to x, so 
Akjh{cc {Qextj — Qext,jr)/a) is independent of the grain radius (except for Fe dust, which 
has Qahs oc a^). For larger grains up to a: < 1, the scattering with Qsca. cx x'^ is more 
important, so Akjh increases. Because, in general, Qext increases to a value of Qext ~ 3 — 5 

= n + ik) is the complex refractive index, and then converges 
iches a maximum near |m — l|a: ~ 2 and then rapidly drops to 


near 

m — 1 

x ~ 2, where 

to 2 ( 

Draine 

2011 

), 


10^ cm^ 


zero. For grains of a = 0.1 pm, Akjh = lO^-lO'^ cm^ g“^. Fe grains have Akjh ~ 
g“^ almost independent of their sizes. As discussed in the previous section, the dust grains 
responsible for the NIR extinction should have large Akjh and/or small FIR Kabs,v compared 
to the general interstellar dust composition. In this sense, small Fe (or other metallic) dust 
grains seem to be a candidate for the cool dust. 


4-2.2. Spectral Modeling 


We now perform spectral htting to the SN dust SED at position A. We assume that 
the SN dust is composed of two temperature components, i.e., warm and cool dust, as noted 
in Section 14.11 A single component cannot explain both the MIR peak and the FIR excess. 
We found that essentially all of the MIR emission is from the warm dust. Therefore, we 
simply assume that the 12 and 24 pm brightnesses are from the warm dust, which then 
yields the temperature and surface density of warm dust without htting the entire SED. 
Color correction was made using the hlter response curve at each band and the spectral 
shape of each grain model. 


Figure [Hl(a) shows the derived temperature and surface density of warm dust for different 
species. Because the FIR nahs,u values of dust grains, except for Fe dust, do not depend 
on their sizes (Figure [6]), the derived temperature and column density of dust species are 
independent of the grain size. Not all of these dust species can match the FIR part of the 
SED. This is shown in Figures[8](c) and (d), where we compare the expected FIR brightnesses 
to the observed brightnesses (the red dashed and dotted lines). Note that it is acceptable 
if the expected brightness is less than the observed brightness because the difference can be 
contributed by cool dust, but not vice versa. Therefore, only three dust species can match 
the FIR part of the SED: MgSiOs, Mg 2 Si 04 , and Si02. The temperature and column density 
of the three grain models are 80 K < T^arm < 120 K and 10" 


g cm ^ < S 


d,warm 


< 10 


-6 


g cm respectively. The compositions and corresponding temperatures of this warm dust 


component are well consi stent with previous observational and theoretical results flRho et al. 


20081: iNozawa et al.ll2010l) . Figure[8](b) shows the expected E{J — H) of warm dust estimated 


from their surface density using Equation [2l Note that the E{J — H) values of the three 
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dust species are less than ~ 10 ^ mag, which is a few orders of magnitude smaller than the 
observed E{J — H), as we already pointed out in Section 1411 

With the warm dust component constrained, we can now search for possible candidates 
for the cool dust component. We perform a two-temperature modihed blackbody £t of the 
SED for each of the three possible warm dust species above, the temperature and surface 
density of which are hxed. The left column of Figure [9] shows the derived temperature and 
surface density of the cool dust for each warm dust species. All combinations can £t the 
observed SED within the la uncertainty of the measured brightness except when the cool 
dust is MgO grains. The opacity of MgO has a small bump peaking at 100 pm, so it cannot 
properly £t the observed brightness at either 100 or 160 pm. The temperature of the cool 
dust ranges from 30 to 70 K depending on its composition. When the warm dust is Mg2Si04, 
the cool dust temperature is systematically lower. The right column of Figure [9] shows the 
expected color excess of each cool dust species. We see that small (< 0.01 pm) Fe grains can 
match the observed E{J — H) 0.2 mag for any warm dust species with Ed = 2 — 5 x 10“^ g 
cm“^. None the other grain species can produce sufhcient extinction. 

Figure IHl however, shows that as the grain size increases, E{J — H) increases because 
larger grains have higher Akjh (Figure |6]), so perhaps grains larger than 0.1 pm might also 
produce E{J — H) ^ 0.2 mag. Figure [7] shows Akjh for each grain species as a function of 
grain radius. For oxide grains (i.e., MgO, MgSiOa, Mg2Si04, AI2O3, and Si02), Akjh reaches 
a maximum at ~ 0.4 pm and then drops abruptly at > 0.6 pm, whereas for metallic grains 
(i.e.. Si, Fe, FeS, and Fe304) and amorphous carbon, the transition occurs at somewhat 
shorter wavelengths. We find that pure Si grains with a radius of 0.16-0.21 pm can also 
explain the color excess we observed because of their large Akjh (Figure [7]), although the 
permitted size range appears to be too narrow. Table [3] summarizes the parameters of the 
warm and cold dust grains that can explain both the FIR SED and the NIR extinction. 


5. Discussion 


5.1. Cool SN Dust Responsible for the NIR Extinction 


According to our result, the cool dust that can explain the observed NIR extinction in 
the southern part of Cas A (i.e., in Slit 4) could be either small (< 0.01 pm) Fe or large 
(> 0.1 pm) Si grains. In this section, we discuss these two possibilities. 


In the Slit 4 position, the re verse shock is now encountering pure Fe ejecta in the 
innermost region flKoo et al.l l2013l) . Considering this, cool SN dust mainly consisting of 
Fe grains seems to be reasonable. The small size also appears to be consistent with the 
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theoretical prediction that the average radius of newly formed SN dust grain s in SN Ilb is 


smaller than 0.01 um. bec ause of low gas density in the expanding SN ejecta flKozasa et ah 


20091: iNozawa et alJl2010h . On the other hand, theory also predicts that the gas density of 
the innermost ejecta is too low to produce pure Fe grains when the temperature of the region 
drops to 800 K, at which Fe grains can start to coagulate, so pure Fe dust is not expected 
to form in SNe lib flKozasa et al.ll2009l: iNozawa et ahll2010l) . However, the above conclusion 
assumes a spherically symmetric explosion with homogeneous and stratihed ejecta. If the 
gas density is inhomogeneous and/or the SN explosion is asymmetric, we cannot rule out 
the possibility of pure Fe grain formation in th e innermost SN ejecta, where the gas density 
is relatively high fjSarangi fc Cherchneffll2014l) . Inde ed, the ejecta rn aterial emitting [Fe II] 
lines is dense, with electron densities of n-u 10"^ cm 


10"^ cm ^ f Koo et ah 20131) . 


Another possible grain species responsible for the large E{J — H) k, 0.2 mag is large 
Si grains. This possibility may be supported by previous infrared studies; Spitzer spectral 
mapping observations found strong [Si II] emission together with [O IV] and [S III] lines in the 
interior of the remnant, which arises fro m the unshocked SN ejecta photoionized by UV and 


X-ray emission from the Bright R im fe.g.. lSmith et al.ll2009l:llsensee et al.ll2010l:lArendt et ah 


—11-—- nhn ' - 

2014j:lMilisavhevic fc Fesenll2015l b The presence of ionized silicon in the interior supports the 
suggestion that the cool dust could be Si and/or silicate grains. Theoretically, pure Si grains 


can f orm in the Si-rich layer, but the grain size is very small, i.e., ;< 0.01 /mi (INozawa et ah 


2010l) . Howev er, again if the SN ejecta is clumped, grains larger than 0.1 /mi can form in 
dense clumps fISarangi fc Cherchneffll2014l) . so large Si grains are not implausible. 


5.2. Composition of SN Dust in Cas A 


From previous studies, it is reasonably well established that there are two SN dust 
components in Cas A: warm (~ 100 K) dust swept up by the reverse shock and cool (30-40 
K), unshocked dust residing in the interior of the remnant. The composition of the warm 
dust is relatively well specihed from its MIR spectral features: Mg-silicates, Si02, AI2O3, 
and/or C glass. The estimated total mass of warm dust ranges from 0.008 to 0.054 Mq. 
The composition of the cool dust has been essentially unknown because in the FIR, the 
SED is mostly smooth without prominent spectral features. Ass uming the FI R absorption 


coeffi cients of the general interstellar dust fjSibthorpe et al.ll2010h or silicates flBarlow et ah 
20101 ). a total mass of ~ 0.07 Mq has been obtained. 


Our analysis also showed that there are warm and cool dust components in the southern 
ejecta shell. The derived temperature (80-120 K) and composition (MgSiOs, Mg2Si04, or 
Si02) of the warm dust are consistent with previous results. For the cool dust, we were able 
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to show for the hrst time that its composition is either small Fe or large Si grains. Their 
temperatures are 30-50 K. These grains are not expected in uniformly expanding stratihed 
SN ejecta, but they can possibly form in dense clumps. 

At hrst glance, the different compositions of warm and cool dust appear to be incom¬ 
patible if we consider that the cool dust is swept up by the reverse shock and subsequently 
turns into warm dust with a temperature of ~ 100 K. However, this can be understood if we 
consider that the SN ejecta is clumped. As we pointed out above, the cool dust responsible 
for the NIR extinction resides in dense clumps. When the dense clumps are swept up by the 
reverse shock, the shock speed is slow (e.g., < 10 km s~^), and the gas temperature is below 
10^ K flKoo et al.ll2013l) . so the dust grains in the clumps are not heated to high tempera¬ 
ture. On the other hand, when the diffuse ejecta is swept up by the reverse shock, the gas 
is heated to X-ray-emitting temperature (~ 10^ K), and the dust grains there are heated by 
collisions with electrons to high temperature. Therefore, it is the shock-heated warm grains 
in the diffuse ejecta that emit most of the MIR and FIR (< 70 /rm) emission (Figure |5]), 
and their compositions could be different from those of grains in dense clumps. In summary, 
our observation can be explained if the SN ejecta is clumped and the grain species in diffuse 
material and dense clumps are different, i.e., silicate grains in the diffuse material and small 
Fe or large Si grains in the dense clumps. The cool dust that we detected could reside either 
in unshocked or shocked clumps, and it is these cool grains that are responsible for the NIR 
extinction. 


When the cool dust is pure Fe grains, we can make a crude estimation of the dust 
formation efficiency from our observation. The required Fe dust mass density of (2-5) xlO”"^ 


g cm ^ corresponds to an Fe nucleus column density of (2-5) xlO 


cm 


For comparison. 


the characteristic size of [Fe II]-line-em itting cluinps is 5" or 0.08 pc at 3.4 kpc, and their 
average electron density is 2 x 10^ cm“^ (IKoo et al.ll2013l) . For pure Fe ejecta clumps, because 
Fe is mostly in Fe"*", it means that the number density of Fe atoms is 2 x 10^ cm“^, and the 
characteristic Fe column density in the gas phase of Fe ejecta clumps is iVpe ~ 5 x 10^^ cm“^. 
Therefore, if there are Udump such clumps along the sight line where E{J — H) ^ 0.23 in Slit 

-3 


4, we obtain a dust-to-gas ratio of 0.4 — 1 x 10 n, 


clump * 


_ ^The total flux of cool dust at 100 /mi in the global SED is ~ 30 ± 10 Jy fiBarlow et ah 

20101 ). The mass has been estimated to be ~ 0.0 7 Mq if we adopt the FIR absorption 


coefficients for silicates from iDorschner et af.l (119951) . The contribution of small Fe or large 
Si grains to this mass budget would depend on the volume filling factor of the clumps but is 
probably not large. To address this issue, one may carry out an analysis similar to that in 
this paper toward the entire ejecta shell. 
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6. Summary 


We performed NIR spectroscopic observations of Cas A in which we obtained the spectral 
and kinematical properties of 63 [Fe II]-hne-emitting knots spread over the main ejecta shell. 
All of the knots show strong [Fe II] 1.26 and 1.64 pm lines, the ratio of which provides 
a direct measure of the extinction. From an analysis of the extinction toward individual 
knots, we showed that the NIR extinction is due in part to SN dust within the remnant. 
We explored the nature of the SN dust responsible for the NIR extinction by analyzing its 
thermal infrared emission. Our main results are summarized below. 

1. We found that [Fe II] emission from redshifted SN ejecta is in general more heavily 
obscured than that from blueshifted SN ejecta (Figure [2]). We interpret the correlation as 
evidence for newly formed SN dust within the remnant. The amount of excess extinction 
varies considerably from one sight line to another, which suggests a highly non-uniform 
distribution of SN dust. 

2. One slit (Slit 4) is located across the southern [Fe II] ejecta hlament, which has the 
velocity structure of an expanding shell. Along this particular sight line, we measured an 
excess NIR extinction of E{J — H) = 0.23 ± 0.05 mag between the front and back sides of 
the shell, which must be entirely due to SN dust (Figure E]). 

3. We analyzed the SED of thermal dust emission toward Slit 4 to show that there are warm 
(~ 100 K) and cool (~ 40 K) SN dust components along the sight line (Figure |5]). Among 
the grain species predicted from theoretical dust formation models in CCSNe, only silicate 
grains (i.e., MgSiOs, Mg 2 Si 04 , Si02) can match the SED of warm SN dust which has a sharp 
24 pm peak. The NIR extinction due to these warm silicate dust grains is negligible, so the 
observed excess NIR extinction should be due to the cool SN dust. 


4. According to our analysis, only two grain species can explain both the SEDs of cool dust 
and the NIR excess extinction: (1) small (< 0.01 pm) Fe grains or (2) large (> 0.1 pm) Si 
grains. However, neither would be consist ent with the dust fo rmation theory for SN Ilb with 
uniformly expanding, layered ejecta fe.g.. iNozawa et al.ll2010l) . which predicts the formation 


of only small (< 0.01 pm) Si grains and no Fe grains. However, if the SN ejecta is clumpy, 
we might expect that Fe grains and large Si grains can be produced in dense clumps. 


5. We suggest that the unshocked SN ejecta is clumpy and the grain species in diffuse 
material and dense clumps are different, i.e., silicate grains in diffuse material and small Fe 
or large Si grains in dense clumps. In the shocked ejecta, both still exist, but it is the shock- 
heated warm grains in the diffuse ejecta that emit most of the MIR and FIR (< 70 pm) 
emission. The cool dust that we detected could reside in either unshocked or shocked clumps, 
and it is these cool grains that are responsible for the NIR extinction. The contribution of 
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each grain species to the total mass (> 0.1 Mq) of cool SN dust in Cas A remains to be 
explored. 
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Fig. 1.— (a) Slit positions (blue bars) marked on the continuum-subtracted [Fe II] 1.64 /rm 
narrow-band i mage obtained in 2008 A ugust. The black cross mark corresponds to the center 
of expansion flThorstensen et al.ll200lh . The green contours represent the hydrogen column 
densi ty (iVn) derived from a single-component £t of Chandra X-ray data flHwang fc Laming 
2012), and the levels are 1.5, 1.9, 2.3, 2.7, and 3.1 x 10^^ cm“^. The red solid line is the 
boundary of the iVn map. (b) Variation of E{J — H) with position angle. The position angle 
is the angle along the perimeter of the remnant measured from north to east with PA = 0 
corresponding to a line drawn due north from the expansion center (black cross ma rk) of 
Cas A. The filled symbols represent the 63 infrared knots detected flKoo et al.l 120131) . with 
green squares indicating He-rich knots, red circles in dicating S-rich kno t s, and blue diamonds 
indicating Fe-rich knots. The seven optical knots of lHnrford fc Feseru fjl996h are marked by 
open symbols: hve FMKs and two QSFs appear a s red circles and green squares, respectively. 
The 19 infrared knots observed in lEriksen et al.l (120091) are also marked by black X symbols. 
The gray-shaded regions are the areas uncovered in our observation. The black solid line 
shows the variation of the average EU — H) along the iron-bright rim obtained from the 
Nn map of lHwang fc Lamingl (120121) . and the black dashed-dotted lines represent maximum 
and minimum E{J — H) at each position angle (see the text for details). 
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Fig. 2.— (a) Variation of — for the 70 knots as a function of radial velocity. Symbols 
and colors of the knots are the same as in Figure [H and the correlation coefficient of the 
knots (p) is given in the bottom right, (b ) Same as (a), bu t after the subtraction of the 
foreground extinction using the map of lHwang fc Lam in d (120121) . 
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Fig. 3.— Position-velocity diagram at (a) 1.26 fim and (b) 1.64 /rm in Slit 4, showing the 
brightness of [Fe II] 1.26 and 1.64 /im emission. The y-axis represents the radial distance 
in arcminutes from the expansion center (see Figure [I]- (a)). The contour levels correspond 
to 2cr and Aa for [Fe II] 1.26 /rm and 3cr and 6cr for [Fe II] 1.64 /mi above the background 
rms noise. Unit of the color bar on the right side of the panels is 10“^^ erg cm“^ s“^ A“^. 
(c) C olor excess (E(J — H)) may of Slit 4. The foreground extinction derived from the iVn 


map flHwang fc Lamindl2m211 is subtracted from the observed E{J — H). The gray hatched 
areas are the areas where the OH line contamination is signihcant compared to the [Fe II] 
emission, (d) One-dimensional E{J — H) prohle of the pixels above 4(T and 6cr rms noise for 
[Fe II] 1.26 and 1.64 /im, respectively. The black solid line represents the intensity-weighted 
extinction prohle, whereas the red, green, and blue solid lines are the mean prohles obtained 
from the slices in (c). The mean velocity resolution and extinction error are indicated by 
the black solid box at the lower right. 
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Fig. 4.— (a) One-dimensional profiles of 12, 24, 70, and 100 pm brightnesses along the length 
of Slit 4. Each prohle is background-subtracted and normalized by the peak brightness at 
lf67. (b) Location of the background region (red box; 10" x 50" between 2f0 and 2'.2 in 
radial distance) on the synchrotron-subtracted PACS 100 pm image. Slit positions are also 
marked with blue boxes. The bottom image is a zoomed-in image of the area in the green 
box in the top frame, and the black scale bar in the lower right represents an angular scale 
of 30". 
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Fig. 5.— (Left) Spectral energy distribution (SED) of dust emission in Slit 4. The red 
and blue solid lines represent the SEDs of total dust continuum emission in regions A and 
B in Figure m and the black solid line is that of the background region (see Table [1]). The 
red and blue dotted lines represent the background-subtracted SEDs in regions A and B, 
respectively, and therefore the SEDs of pure SN dust. The symbol with a downward arrow 
denotes the la upper limit of the brightness. (Right) Two-component ht for the SED at 
position A. The warm dust (red dashed line) is assumed to be MgSiOs, whereas the cool dust 
(blue dashed line) is assumed to have the general interstellar dust composition (see the text 
for details). Red and black symbols represent the measured and color-corrected brightnesses, 
respectively. 
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Fig. 6.— i^abs,u at 24 /xm (Kabs, 24 ; upper panel) and 70 /im (Kabs, 7 o; lower panel) versus 
Akjh for 10 grain species in Table [2l Symbols and colors denote the grain species and sizes, 
respectively. We connec t the same sy mbols with black dotted lines. For comparison, the 
general interstellar dust (iDrainel 120031) is marked with a star symbol. 























29 



Grain Radius (urn) 


Fig. 7. 


^i^jH versus grain radius for 10 grain species. 
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Fig. 8.— (a) Temperature and column density of warm dust grains obtained from the fit of 
MIR emission. Symbols and colors are the same as in Figure El (b) Predicted E{J — H) of 
warm dust grains from their column densities in (a). The red dashed and dotted lines are 
the measured color excess and its la uncertainty, respectively, (c) Expected brightnesses of 
warm dust grains at 70 and 100 pm. The red dashed and dotted lines represent the observed 
brightnesses and their la uncertainty, respectively, (d) Same as (c) but for 160 and 250 pm 
brightnesses. 
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Fig. 9.— (a) Temperature and column density of various cool dust species obtained from 
the fit of FIR emission assuming that the warm dust component is MgSiOs. (b) Predicted 
E{J — H) of cool dust from the column density in (a). The red dashed and dotted lines are 
the measured color excess and its la uncertainty, respectively, (c)-(d) Same as (a)-(b) but 
assuming the warm dust component is Mg 2 Si 04 . (e)-(f) Same as (a)-(b) but assuming the 
warm dust component is Si 02 . 
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Table 1. Observed and dust-continuum infrared brightnesses in Slit 4 


Waveband 


Observed 

a 


Dust-continuum 


A 

B 

Background 

A 

B 

Background 

WISE 12 

24.0±1.1 

7.9±0.4 

6.6±0.5 

21.9±1.0 

7.5±0.3 

6.5±0.5 

MIPS 24 

397±16 

37.1±1.5 

34.6±5.3 

378±15 

31.0±1.2 

33.9±5.3 

PACS 70 

170±8 

74.5±3.7 

68.8±5.7 

163±8 

74.5±3.7 

68.8±5.7 

PACS 100 

210±11 

175±9 

151±11 

201±10 

175±9 

151±11 

PACS 160 

296±15 

285±14 

266±18 

296±15 

285±14 

266±18 

SPIRE 250 

205±14 

206±14 

200±15 

205±14 

206±14 

200±15 

SPIRE 350 

115±8 

113±8 

113±8 

115±8 

113±8 

113±8 

SPIRE 500 

48.0±3.4 

47.0±3.3 

48.5±3.4 

48.0±3.4 

47.0±3.3 

48.5±3.4 


^The observed brightnesses include contribution from emission lines. 


Note. — Brightnesses in units of MJy sr The uncertainties are derived from a quadra¬ 
ture sum of the photometric uncertainty and the brightness variation within the aperture. 
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Table 2. Dust grain species considered in this work 


Species 

Condition^ 

Type 

Density 
(g cm“^) 

A Coverage References 
(^m) 

C 

u 

Amorphous 

2.28 

4E-2 

- 2E3 

1 

MgO 

u 

Crystalline 

3.59 

2E-3 

- 625 

2 

MgSiOa 

m/u 

Glassy 

3.20 

2E-1 

- 500 

3 

Mg 2 Si 04 

m/u 


3.23 

lE-1 

- 1E5 

4 

AI 2 O 3 

m/u 

Amorphous 

4.01 

2E-1 

- 500 

5 

Si 

u 

Amorphous 

2.34 

lE-2 

- 148*’ 

6 

Si02 

m/u 

Amorphous 

2.66 

lE-4 

- 500 

7 

Fe 

u 


7.95 

lE-1 

- 1E5 

4 

FeS 

u 


4.87 

lE-1 

- 1E5 

4 

Fe304 

m 


5.25 

lE-1 

- 1E3 

8 


“‘u = unmixed SN explosion model, m = mixed SN explosion model 

'^The maximum of A coverage is limited to 148 /im. In this case, we assume 
that the optical properties, n and fc, do not vary at longer wavelengths. 


References. — ( 11 IZubko et al.l (Il996ll. (2 ) iRoessler fc Huffman 


(3) Dorschner et al. ( 1995h. Semenov et al.1 (2003), (5) 


< 8 


/xm 


for A 

(Il985h . (7) iPhiliPD 


and 


Koike et al 


(119911) 

(I995I1 


Beeemann et al.l (|1997l) for A > 8 fxm, ( 6 ) IPiller 


(|l985l) . ( 8 ) A. Triaud and H. Mutschke (unpublished; see 


http://www.astro.uni-jena.de/Laboratory/DCDB/mgfeoxides.html). 
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Table 3. Parameters of SN dust in Slit 4 




Warm Dust 



Cool Dust 


Species 

T 

Sd 

E{J - H) 

Species 

T 

Sd 



(K) 

(10“® g cm“^) 

(10“^ mag) 


(K) 

(10“^ g cm“^) 

(mag) 

MgSiOa 

(a = 0.1/im) 

115 ±2 

0.032 ± 0.005 

0.11 ±0.02 

Si 

(0.16 < a < 0.21/.tm) 
Fe 

(a < 0.01/xm) 

45 ±6 

0.039 ± 0.022 

0.15 - 0.65 




< 45 ± 7 

> 2.0 ± 1.3 

> 0.16 ±0.10 

Mg2Si04 
(a = 0.1/im) 

92 ± 1 

0.14 ±0.02 

0.60 ±0.09 

Si 

(0.16 < a < 0.21/im) 
Fe 

(a < O.Ol^m) 

30 ±7 

< 26±6 

0.054 ±0.052 

> 6.9 ± 7.1 

0.21 - 0.90 

> 0.55 ±0.56 

Si02 

(a = 0.1/im) 

78 ± 1 

0.9 ± 0.1 

2.2 ±0.3 

Si 

(0.16 < a < 0.21/im) 
Fe 

(a < 0.01/xm) 

42 ±7 

< 40±7 

0.035 ±0.025 

> 2.1 ± 1.8 

0.14 - 0.59 

> 0.17 ±0.14 









